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Abstract 

00 

CN ' We obtain the ANNNI model from a Heisenberg model with large single-ion 



anisotropy energy, D, as might be relevant for helical spin systems. We treat 

quantum fluctuations to lowest order in 1/5 at zero temperature within an 

expansion in J/D, where J is an exchange energy. The transition from the 

state with periodicity p = 4 to the uniform state (p = oo) occurs via an 

infinite sequence of first order transitions in which p increases monotonically. 
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Systems with long-period modulated structures are surprisingly common in nature. Ex- 
amples include helical phases in the rare-earths and their compounds [ [I], polytypismf Q, 
and the arrangement of antiphase boundaries in binary alloysf |c|. A given compound may 
exhibit many different modulated structures of differing wavelength as a control parameter 
such as the temperature is varied. Some modulated structures can usefully be viewed as an 
assembly of domain walls when the energy for introducing a wall passes through zero. The 
stability of the different structures is then determined by the interactions between pairs, 
trios, etc. of walls [ §]]. It has been established that these interactions can result from en- 
tropic contributions to the free energy [ H and from softening of the spins [ |J. Here our aim 
is to show that quantum fluctuations can also stabilize long-period modulated structures. 

The Hamiltonian we consider is 
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where i labels the planes of a cubic lattice perpendicular to the z-direction and j the position 
within the plane. Also (jj') indicates a sum over pairs of nearest neighbors in the same 
plane and Sjj is a quantum spin of magnitude S at site (i,j). For D = oo, only the states 
Siz = &iS, where Oi — ±1 are relevant and TC reduces to the axial next-nearest neighbor 
Ising (ANNNI) model, first proposed to describe helical phases of the heavy rare earthsf 0, 

ft A = —Jo 2^ Oij°ij' ~ J\ / iVij&i+ltJ 

+ J2j2 <T iJ <T i+2>3 ■ ( 2 ) 
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The ground state of the ANNNI model is ferromagnetic for k = 3%j J\ < 1/2 and an antiphase 
structure with layers ordering in the sequence {<7j} = {...1,1, — 1, — 1, 1, 1, —1, — 1 . . .} for 
K > 1/2. k = 1/2 is a multiphase point [ [|], where the ground state is infinitely degenerate 
with all possible configurations of ferromagnetic and antiphase orderings having equal en- 
ergy. For classical spins S = oo, the ground state (and therefore the multiphase point) is 



maintained as D is reduced from infinity as long as D is larger than about 1/2. For higher 
order anisotropies this is not the case[ ||. 

To describe how the degeneracy is broken at the multiphase point we use a notation 
similar to that of Fisher and Selkef [| so that (m, ri2, . . . n m ) denotes a state in which spins 
form domains (of parallel spins) whose widths repeat periodically the sequence n\, n 2 , . . . n m . 

Fisher and Selke[ || showed that at nonzero temperature T the degeneracy at the mul- 
tiphase point is broken to give a sequence of phases (2 fc 3), for k = 1,2,3. . .. Fisher and 
Szpilka[|]|| later recast their analysis in terms of domain wall interactions and we will follow 
their formulation. 

In view of this interesting phase diagram in the k-T plane, we are led to study the 
phase diagram in the k-D^ 1 plane when the spins are quantum operators. That quantum 
fluctuations can remove ground-state degeneracies was pointed out by Shender[|9| and given 
the apt name "ground state selection" by Henley [ [I0| . In this paper we show how the 
multiphase degeneracy is resolved by quantum fluctuations. 

To study quantum fluctuations we introduce the Dyson-Maleev[ [TTJ] transformation 
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where 5 a ^ is unity if a = b and is zero otherwise and af (oj) creates (destroys) a spin 
excitation at site i. We thereby transform the Hamiltonian of Eq. (Q) into the bosonic form 



Hda,}) =E + H + Vn + V v + V& 
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where E = Ha, 
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with D = D + 2 Jo and Vu (Vw) is the interactions between spins which are parallel (antipar- 
allel) 
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where X(i, i';j) \Y(i,i r ;j)] is unity if spins (i,j) and (i',j) are parallel [antiparallel] and is 
zero otherwise. In Eq. (|4]) V^' represents the four operator terms proportional to 1/S 2 . 
Fluctuations out of the classical ground state (the boson vacuum) only occur at the walls due 
to Vj|/. We do not consider quantum fluctuations within a plane, since the phase diagram 
is determined by the interplanar quantum couplings. Also, since the walls in this three 
dimensional system are flat at T = 0, we may characterize states of the system in terms of 
distances between walls. 

We now consider the structure of perturbation theory for all states which are degenerate 
at the multiphase point k = 1/2. Perturbation theory generates corrections to the diagonal 
energy of the classical states in powers of 1/S and J/D, where J = J\ or J 2 . Off-diagonal 
matrix elements (for example, in which two domain walls both move through one lattice 
constant) first occur in 25'th order perturbation theory and may be ignored. We will only 
include effects of the quadratic Hamiltonian, i.e. we will work to leading order in 1/S. 

Instead of a direct evaluation of the energy of all possible phases, we follow the methods 
of Fisher and Szpilkaf |8| and study the sequence of wall interaction energies: E w , the energy 
of an isolated wall; V 2 (n), the interaction energy of two walls separated by n sites; and 
generally Vk(ni,n 2 , . . .rik-i), the interaction energy of k walls with successive separations 
rii, n 2 , ... Uk-i- In terms of these quantities one may write the total energy of the system 
when there are walls at positions rrii as 



E = Eo + n w E w + J2V 2 {mi +1 -mi) 

i 

+ Yl Vs(™*+2 - m i+1 ,m i+1 - rrii) 

+ Y Vi(m i+3 - m i+2 , m i+2 - m i+1 , m i+1 - mi) 

+ -'- , (8) 

where E is the energy with no walls present and n w is the number of walls. The scheme of 
Ref. [ [8] for calculating the general wall potentials Vu is illustrated in Fig. 1. Let all spins 
to the left of the first wall have <Ji = o and those to the right of the last wall have <Ji = 77 for 
k even and oi = —77 for k odd. The energy of such a configuration is denoted E^a, 77). If 
a = —1 (77 = —1) the left (right) wall is absent. Then the energy ascribed to the existence 
of k walls is given by[ 

Vk(ni,n 2 ,...n k -i)= Y ar)E k ((7,r]) . (9) 

a,r)=±l 

Contributions to E k which are independent of o or 77 do not influence V&. Ej,(<j,r]) is 
calculated by developing the energy in powers of the perturbations V\\ and Vy. To lowest order 
in 1/D, contributions to V& can be obtained, for instance, by creating an excitation at the left 
wall (using Vw) and (for wall separations n\ > 3) using VTi to hop the excitation sufficiently 
near the other wall that one (or more) energy denominator depends on 77. Examples of such 
processes are shown in Fig. 2. 

For instance for the top diagram of Fig. 2, we get 

E 2 (a,rj) = - (— ^ ) %i , (10) 

which gives a contribution to V 2 (2) at order J 3 /D 2 S of 
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Collecting all such processes we find the general result 

V 2 (2n) = ^'WAljM+SJ, ^ (13) 



These results may be understood in terms of a correlation length £ ~ [1 / ln(4D / J 2 )} which 
governs wall-wall interactions. 

More generally, power counting shows that 

V 3 (2n, 2n) ~ V 3 (2n, 2n + 1) ~ J(J/D) 4n 
V 3 (2n - 1, 2n - 1) ~ V 3 (2n - 1, 2n) ~ J(J/D) 4n ~ 1 (14) 

and Vfc(ni, n 2 , . . . n^-i) ~ J(J/D) X , where x > J2j n j ~ 2. Second order perturbation theory 
yields the result 

E w = 2Ji - 4J 2 - ^tM + ( J 3 /^ 2 ^) • (15) 

When E w > 0, the ferromagnetic phase is stable. This happens for J 2 < J c = J\/2 — 
{3Jf/8DS) . . .. 

We wish to describe the sequence of phases which occur as J2/J1 is decreased starting 
from (2) when J2/J1 > 1/2 and reaching (00) when J 2 < Jc- As Fisher and Szpilka show, 
the phase boundary along which (n) and (n + 1) have the same energy is given by 

E w = nV 2 (n) - (n + l)V 2 (n + 1) + nV 3 (n, n) 

-(n + l)V 3 (n+l,n + l) + ... (16) 

This relation yields a critical value of J 2 , denoted J nc which can be expressed as J nc = 
J c + AJ 2 (n), where 



AJ (n) _ nKjn) - {n + l)V 2 (n + l) 
21 ' ~ dE w /d.J 2 
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J 2 =Ji/2 



Thus, to elucidate the topology of the phase diagram, it is not necessary to know J c accu- 
rately. For n not too large, Eqs. ([T2|), ([13|), and (0) give AJ 2 (n) ~ V 2 {n) ~ J(J 2 /D 2 )^ 2 \ 
where [x] is the integer part of x. The tentative conclusion is that one has successive regions 
of stability of the phase (n), where n increases as J 2 decreases, as shown in Fig. 3. However, 
we must check the stability of the phase boundary to mixed phases of (n) and (n+1). 



As Fisher and Szpilka show, the condition that this phase boundary be stable is that 
F n < 0, where 

F n = V 3 (n, n) - 2V 3 (n, n + 1) + V 3 (n + 1, n + 1) . (18) 

Here the last term is higher order in 1/D than the first two and can be neglected. All per- 
turbative terms which contribute at lowest order in 1/D to V 3 (n,n + 1) have their analogs 
for V 3 (n, n). By an appropriate grouping of terms one can show that for n > 2, F n < 0. Ba- 
sically this happens because even order ground-state-to-ground-state terms in perturbation 
theory are negative. The case n = 2 is special in that F2 = at lowest order. Then it is 
necessary to go to the next order, where we find 

V 3 (2,2) = jtfjLl-J? + ^-h - 2J 2 2 ] 

12 J 2 

+ ~ 2 [-4J? + 12J 2 J 2 - 5Ji J 2 2 + 10Jf] (19) 

o t4 19 72 

\/ 3 (2, 3) = -d^- + J±!*-[2J*J 2 + 4 Jx J 2 + 5J|] . (20) 

V ' ; (4D) 3 5 (4 J D) 4 ,S L ' 2 2J V ; 

\/ 3 (3,3) = 0(J 6 / J D 5 ,S). (21) 

To leading order in 1/S we may set J 2 = Ji/2, in which case the above results indicate that 
F 2 ~ A/D 4 , where A < 0. Thus all the phase boundaries between phases (n) and (n + 1) 
are stable agains subdivision. 



The above results are valid (as we shall see) for n <^ y D/J. When this limit is violated, 
the entropy of more complicated perturbation contributions can compensate for taking more 
powers of J/D. We overcome this limitation with respect to V 2 (n) as follows. We work to 
lowest (second) order in Vw (A pair of excitations is created, one to the left of the left wall 
and one to the right of the left wall, as in Fig. 2, and is later destroyed.) To simplify the 
result we assume that the excitation created to the left of the left wall does not propagate. 
We work to first order in the field exerted on spins n — 1 and n by the spins in the neighboring 
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domain. The result for the ground state energy is then expressed in terms of the EXACT 
spin- wave Green's function, G^ n \ for an isolated domain of n spins. In this way we sum over 
all trajectories of the spin deviation inside the domain of n parallel spins. For small n we 



reproduce the above resultsf T3[. For large n the result at leading order in J/D is 

V 2 (n) = 4J|G (n) (2,n - if / S , (22) 

where 

G'lul^^. (23) 

Here <p a and e a are the exact eigenstates and energies for the single-spin excitations of an 
isolated system of n parallel spins. We carried out an exact evaluation of G^ n \2, n — 1). For 
large n, we found 

AD 

w 

'2 



G(2,n- l) 2 = ^ye n ^sm 2 (n5) , n even 

^2 

AD 



^-e~ n/ t cos 2 {n5) , n odd , (24) 

where 5 = Ji/yl6DJ 2 and £ — £ differs from zero due to corrections which are higher order 



in J/D. Eqs. (22"D and ([M|) seem to imply that V^n) can become arbitrarily close to zero. 
That is an artifact of truncating these equations at leading order in J/D. Where V^n) 
would be small, one must keep the appropriate terms which are otherwise corrections. So 



doing we have a result which is uniformly asymptotically correct for n ^> J D / J[ |i~4"f : 

, . lbD _„/?/ . 9/ c- i\ J1J2 \ 
\ / 2{ n ) — — tt~ e sm \ n v — I— <a>) H =— , n even 



S V v 2D 2 

^e-^(cos 2 M + 0) + ^),nodd, 

(25) 



where is a phase shift of order 1/vD. 

An elegant graphical interpretation of the phase boundaries suggested by Fisher and 
Szpilka is that one should construct the extremal convex envelope of V^n) versus n. The 
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points [n, V^n)] which make up the envelope correspond to the phases (n) which occur when 
V2(n) is not convex, as Eq. (|2^ ) shows to be the case. As a result, we conclude that there 
is an infinite sequence of phase boundaries. When 2n5/ir is nearly an integer, the phase 
boundaries will be between phases n and n + 2 because of the nonconvexity of V^n). Note 
that in contrast to the ANNNI model, here V^n) does not pass through zero. This difference 
can be understood as follows. In the present model in order for an excitation to sense the 
presence of a second wall, it has to travel from one wall to the other wall and return, giving 
rise to the factor G 2 in Eq. (22). In the ANNNI model the analogous factor involves only 
a one-way connection corresponding to G. As a consequence, for the ANNNI model the 
sequence of phases terminates at a value, n , which diverges as T — > 0. There is no cut-off 
on n in the present model. 

We were unable to carry out a precise analysis for V3 at large n. Accordingly, at large 
n we can not guarantee the stability of these phase boundaries. It is conceivable that our 
result for small n breaks down and that the phase boundaries obtained from V^n) become 
unstable to mixing, which could even be hierarchical. 

To summarize: 1) We have shown that quantum fluctuations do remove the infinite 
degeneracy of the multiphase point of the ANNNI model. 2) We have shown that quantum 
fluctuations at T = lead to a sequence of first order transitions similar to that for the 
ANNNI model, but involving a different sequence of phases. 3) In contrast to the ANNNI 
model there is no cut-off at large n on the appearance of phases because here V^n) never 
becomes negative. As we explained, this is a peculiarly quantum effect. 
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FIGURE CAPTIONS 

FIG. 1 Configurations needed to calculate the interaction energy for two walls at separation 

n (top) and three walls at separations n and m (bottom). 

FIG. 2 Examples of configurations needed to calculate V^(2) (top), 1^(3) (middle), and 

1^(4) (bottom). Here "+" ("-") indicate creation (destruction) of a spin excitation and the 

arrow indicate a hopping using Vm. 

FIG. 3 Schematic phase diagram of the "soft" ANNNI model. The phase boundary between 

(n) and (n + 1) depends on a power of 1/D which increases with n. We did not attempt to 

represent this dependence on D correctly. 
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FIG.3 
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